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Due to the great level of wind energy penetration in the existing network, 
huge efforts have been directed to enhance the grid-connected wind 
generator performance. This paper shows an application of a hybrid 
algorithm of the particle swarm optimization and the gravitational search 
algorithm (PSO-GSA) to enhance the transient stability of the grid-tied wind 
energy conversion system. The variable-speed wind turbine (VSWT) direct- 
drive permanent-magnet synchronous generator is connected to the network 
through a full-scale converter. The generator- and grid-side converters are 
controlled by utilizing an optimum proportional-integral (PI) controller. The 
criterion of the integral squared error is utilized as an objective function. The 
PSO-GSA based-PI controller efficacy is validated by comparing its results 
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with that are obtained by utilizing the genetic algorithm (GA)-based-PI 
controller. The performance of the suggested control scheme is checked 
during various fault conditions. The control scheme quality is legalized by 
the simulation results that are obtained using MATLAB/Simulink program 
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1. INTRODUCTION 

Recently, massive attention has been directed to renewable energy around the globe, for various 
reasons, such as the cost growth of the traditional power generation, environmental needs, and probability of 
the demise of fuel. According to the latest statistics, it is expected that new installations of more than 55 GW 
will be added in the grid each year until 2023 [1]. Based on the enormous level of wind energy participation 
in the existing electrical grid, some problems have arisen. Some of these problems are the fluctuations of the 
wind farms out put power and the low voltage through capability. So, the wind energy transient 
characteristics have to be studied. 

The variable-speed wind generator system (VSWG) has many merits than that of the fixed in the 
wind energy conversion system (WECS). Some of these merits are the capability to capture maximum power 
and reducing the mechanical stress on the drive-train [2, 3]. There are various types of electrical generator 
used in the VSWG, such as switched reluctance generator (SRG) [4], a doubly-fed induction generator 
(DFIG) [5, 6], and permanent magnet synchronous generator (PMSG) [7-13]. The PMSGs have many 
advantages, for example, they have high number of poles, therefore there are power can be achieved at low 
rotational speed and there is no needing for gearbox resulting in reducing cost. 

The grid-tied variable-speed wind turbine (VSWT)-PMSG is composed of both full-scale frequency 
converters at the generator- and grid-side, intermediate capacitor (C), transformer as step-up, and two 
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transmission lines. The vector control has utilized in controlling the generator- and grid-side converters. In 
this study, the vector control depends on the proportional-integral (PI) controller [14-17]. The proportional- 
derivative (PD) and proportional-integral-derivative (PID) controllers are not more attractive to be utilized in 
the industry due to the disadvantages of derivative control action, where it increase the input harmonics 
frequency to the system, and to prevent this problem, a designed filter should be utilized. Therefore, PI 
controller is commonly used in the industry. PI controllers have characterized by having a large stability 
limit, although they are affected by the variation of the parameters. Four PI are utilized for each converter, so 
the fine-tuning of these controllers is arduous, especially in the non-linear system. 

Many approaches are utilized to reach the optimal parameters of the PI to increase the efficiency of 
the grid-VSWG. Traditionally, Newton-Raphson (NR) and Ziegler-Nichols (ZN) approaches are used to fine- 
tune the PI controller, but these optimizations methods are affected by the type of solver and the initial 
conditions. Besides, the heuristic approaches, such as the gravitational search algorithm (GSA), a genetic 
algorithm (GA), the particle swarm algorithm (PSO) [18], a grey wolf optimizer (GWO), and the ant colony 
optimization (ACO) are used to obtain the optimal values of the PI controllers. But these methods suffer from 
some drawbacks in finding the ideal solution. The GA technique can be applied for both continuous and 
discrete parameters, but it need a complex analysis, a long time, and can be trapped into local optima. The 
PSO has some advantages such as its simplisty and it has few parameters to adjust. But, PSO can have 
difficulty to get the initial design parameters, and it can converge early and can be trapped into local 
optimum solutions with complex problems. GSA depends on Newton’s law of mass attraction and gravity. 
GSA has some merits, including the ability to solve highly nonlinear optimization problems and has stable 
convergence characteristics. But, GSA has drawbacks such as it has slow searching speed in the last 
iterations and has long computational time. ACO can give rapid discovery of reasonable solutions and has 
guaranteed convergence. But, ACO has an uncertain time to convergence and has a complicated theoretical 
analysis. The GWO has some demerits, such as need to long time in convergence and inadequate local 
searching ability. 

The main characteristics of the optimization process are the balance between the exploration and the 
exploitation to reach to the optimal solution. This article presents the application of PSO-GSA [19-23] based- 
PI to control the frequency converter (FC). The PSO-GSA has many merits, such as better capability to 
escape from local optima solution and has a faster convergence rate. The performance of PSO-GSA is 
ensured by comparing the results are achieved by PSO-GSA based-PI controller with that obtained from the 
application of a GA-based PI controller. The various types of faults are considered in this comparison. These 
simulation results are performed by utilizing the MATLAB/Simulink programming. 


2. MODEL SYSTEM 

The main component to adjust output of the VSWT-PMSG is the FC, as shown in Figure 1. The FC 
is composed of two full-scale power converters at the generator and grid-side, C. the capacitor is protected by 
an over-voltage protection scheme (OVPS) during the fault. Each vector control scheme has four PI 
controllers to control the generator and grid-side converters. The generator-side converter (GSC) is managed 
to obtain maximum power at power factor equal one, from the generator. The terminal voltage at the 
common coupling point (Vpcc) and the DC voltage (Vpc) are regulated by the grid-side inverter (GSI). The 
interconnection between the output of the VSWT-PMSG is composed of transformer and double circuit 
transmission line. The data of PMSG are illustrated in Table 1. 
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Figure 1. System model 


Table 1. Generator parameters 


Prasad V,ateq Frequency No. poles 
15MW 575V 20Hz 48 

Rstator Xa Xq Field Flux 
0.006 pu 0.7pu 0.3 pu 0.9 pu 
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3. WIND TURBINE MODEL 
The relation among the speed of the wind and mechanical power can be presented as follow [3]: 


P, = 0.5 prR? V} Cp(à, B) (1) 
where V, is the speed of the wind, p is the density of air, R is the radius of the blade, Cp presents the power 


coefficient, which depends on turbine pitch angle, f and tip speed ratio, A and where w, means the rotational 
rotor speed it can be modeled as: 


Uae Q) 
Leo 5 

Goap) (Grea) 3) 
Cp (å, B) =0. 73( (=) — 0.588 — 0.002(B214) — 13.2) e7184/Li ids 


VSWT is preferred over the FSWT in WECS due to it able to amplify the extracted power at the wind 
speed variations. Figure 2 exhibits turbine characteristic. It is noted that there are different curves for the 
extracted power at various wind speeds. For each curve, a specific point where the maximum power can be 
captured. So, the adjustable controller has to follow the particular rotational speed to get the maximum power 
at each wind speed. This scheme is known as maximum power point tracking (MPPT) 

To get Pnax in terms of w, utilize this equation [3]: 


Or 


R 
Prax = 0.5p7R?( )Cp-opt (5) 


on 


where Cp-opt is coefficient of optimum power and A,,; is optimum value of the tip speed ratio. The rated 
power will be otherwise called P,,,, which is used as the adjusting point for the power at FC. 
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Figure 2. Wind turbine characteristic with MPPT curve 


4. THE FC CONTROL 

The controlling the behavior of the grid-connected VSWT-PMSG is based on the FC, to extract 
Pmax from the generator and adjusting theVpcc, Vpc at various operating conditions. It composes of the GSI, 
GSC, and C. 


4.1. GSC 

The GSC has the major rule to deliver Pnay from the generator at power factor equal one. This 
output is obtained by adjusting the dq -axes components of the stator current (Ig, Iq). The Pmax from the 
generator is achieved by adjusting the J, according to the set-point (Popt), which is obtained from MPPT. 
The J, is adapted to deliver the power when the power factor equal one. As illustrated in Figure 3, there are 
four PI-controllers to control the GSC to acheive the dq-axes components of the reference voltage (V" qq). 
These components are transformed into three-phase components (V* a,b,c). A comparison is occurred between 
V* ab, and a 1 kHz triangular carrier waveform (TCW), to get firing pulses (FP) for the converter. 
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six pulses to 
IGBT switches 


Carrier signal 


Figure 3. The generator-side converter control block diagram. 


4.2. GSI 

The GSI has the responsibility to control Vpç and Vpcc by adjusting the d-axis and q-axis grid 
current components (Ign .Jgn). The Ign is adjusted to control of the Vpc, and the Ian is adjusted to maintain the 
Vpcc to be stable at sever disturbances. In this analysis, the C is selected as 10 mF with a rated voltage with 
rated value equal to 1.15 kV. As obtaned in Figure 4, the grid side inverter control block diagram, to extract 
the reference voltage (V*an,qn) by utilizing four PI-controllers. The V*4n,qn are transformed into three-phase 
signals V* an bn,cn by utilizing a transformation angle (0+). A comparison is occurred between V* an,bn,cn With 
a TCW with 1050 Hz to produce FP to the GSI. To extract (0+), a phase-locked-loop (PLL) is utilized. 


six pulses to 
IGBT switches 


Carrier signal 


Figure 4. The grid-side inverter control block diagram 


4.3. OVPS 

During the fault condition, the reactive power (Q,.) is delivered to the grid. At the same time, no real 
power (P,) is transmitted to the network, which results in increasing the Vpç and leading to destroy the C. The 
OVPS is needed to protect the C from damaging. The components of the OVPS are the controlled IGBT with 
series in resistance [24]. When Vpc exceeds 20% above the nominal value of the Vpc, IGBT will operate to 
dissipate the power through the resistance. So, OVPS will prevent the damaging of the intermediate C. 
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5. HYBRID PSO-GSA AND PROBLEM FORMULATION 

This optimization method is the merging of two algorithms GSA and PSO. The combination of the 
GSA and PSO is performed by utilizing the low-level co-evolutionary heterogeneous hybrid. PSO-GSA is 
defined as co-evolutionary as it runs the two algorithms in parallel combination. The development of the 
PSO-GSA depends on the capability of the PSO for Social thinking and GSA ability to get the local search. 
The merging can be mathematically modeled as follow: 


Vit + 1) = w * V(t) + c; * rand * ac; (t) + c2(t) * rand * (Gpest — Xi) (6) 
where V;(t) shows the i agent velocity at t iteration, ac;(t) is the agent i acceleration at iteration t, cj is a 


weighting factor, w is a weighting function, pest is the best solution so far. At each iteration end, the particle 
velocity and position are updated as follow: 


Xit +1) = X(t) + Vit + 1) (7) 


For the PSO-GSA process, firstly there is initialization for all agents randomly. Then for all agents, 
the resultant force among each other, gravitational constant, and gravitational force are calculated. After this 
phase, the GSA is utilized to get the acceleration of the particles. The best solution has updated at the end of 
each iteration. Then by utilizing the equations (6, 7), for all agents, the velocity and the position can be 
calculated. The updating process can be finished, when the criteria of the velocity and the position reach the 
end. The search space is explored by the agents are founded close to the real best solution; they attract other 
agents. The exploiting of the global best depends on the agents are near to the best solution and move very 
slowly. The cl and c2 are regulated to obtain the equilibrium among the local search and the global search. 
Figure 5 shows the flowchart of GSA. Figure 6 presents the steps of PSO. Figure 7 presents the flowchart of 


PSO-GSA. 
Initialize the parameters of PSO 


( First swarm generation ] 


Evaluate the particles fitness 
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Final global 
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Update the velocity 
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end of 
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Return best 
solution. 


Figure 5. Flow chart of gravitational search Figure 6. Flow chart of particle swarm optimization 
algorithm 
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Figure 7. The hybrid PSO-GSA algorithm flow chart 


Choosing of the objective function in our problem is done by using integral square error (ISE) 
criterion, as expressed: 


ISE = f,((er? +e? + e3” +. e47) + (es? + e6? + €7” + €g”)). dt (8) 


where e4, e2, @3, and e, are the actuating errors Of Pont, lq; Qpmsg» and Ig, respectively, for generator-side 
and es, €e, €7, and eg are the actuating errors of Vpc, Jan, Vecc, and Ign, respectively, for the grid-side. 


6. SIMULATION RESULTS 

The simulation results are achieved by using the MATLAB/Simulink program. The simulation time 
step is 5 micro sec. The results of the PSO-GSA based-PI are compared with that are produced by the GA-PI 
to confirm the PSO-GSA based-PI efficacy under the various disturbances. The three-line to ground fault 
(3LG) is occurred at 1.7 sec at point F, as shown in Figure 1. At 1.8 sec, the circuit breaker at the faulted line 
is opened and then reclosed at 2.1 sec after clearing the fault. Figure 8 (a) represent the response of the P. 
utilizing PSO-GSA-PI controller is faster and better damped than that achieved by GA-PI. Figure 8 (b) 
illustrates the response of Q,. It is noticed that the response of Q, by PSO-GSA based-PI has lower settling 
time than that is achieved by GA-based-PI. Figure 8 (c) illustrates the Vpc¢. Notably, Vpc has lower overshoot 
and better-damped response by utilizing PSO-GSA based-PI than that is obtained by utilizing the GA-PI. 
Figure 8 (d) presents the Vpcc response. Notably, at the steady-state, the response of Vpcc using PSO-GSA 
based-PI has low settling time and has a small error than that obtained by utilizing GA-PI. 
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Figure 8. 3LG temporary fault responses (a) P.. (b) Q,. (c) Voc. (d) Voce 
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For further proof of the PSO-GSA based-PI controller efficiency, the system is subjected to 
unsymmetrical faults. Single line-to-ground fault (1LG), the line-to-line fault (LL), and double line-to-ground 
fault (LLG). Figure 9 shows the responses of the Vpcc according to unsymmetrical faults. Figure 9(a) 
illustrates the Vpc¢ at the LLG fault. Figure 9 (b) shows the Vpcc at the LL fault. Figure 9 (c) shows the Vpc¢¢ 
at the 1LG fault. Finally, the PSO-GSA based-PI can improve the transient stability of VSWT-PMSG under 
various disturbances. 
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Figure 9. Responses of Vpccunder the unsymmetrical faults, (a) LLG, (b) LL, (c) 1LG 


7. CONCLUSION 

This article has illustrated an implementation of the PSO-GSA based-PI to improve the behavior of 
the grid-tied VSWT-PMSG. The PSO-GSA based-PI is legalized by comparing its results with that are 
achieved by utilizing GA based-PI controller. The suggested strategy is utilized in controlling the FC to 
deliver and smooth the P,, and to adjust the terminal voltage to be constant at various operating conditions. It 
is clear from the simulation results that the transient system responses using PSO-GSA based-PI controller 
have better-damped responses and minimum oscillations than that are obtained using GA-based PI. The 
efficiency of the PSO-GSA based PI is checked under various fault conditions. Finally, the PSO-GSA based- 
PI controller has efficiency to deal with controlling the non-linear system. Therefore, the proposed control 
scheme can be used to achieve enhanced response in other systems, smart grids, and storage energy systems. 
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